Quantization of adiabatic pumped charge in the presence of superconducting lead 
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We investigate the parametric electron pumping of a dou- 
ble barrier structure in the presence of a superconducting lead. 
The parametric pumping is facilitated by cyclic variation of 
the barrier heights xi and X2 of the barriers. In the weak cou- 
pling regime, there exists a resonance line in the parameter 
space {xi,X2) so that the energy of the quasi-bound state is 
in line with the incoming Fermi energy. Levinson et al found 
recently that the pumped charge for each pumping cycle is 
quantized with Q — 2e for normal structure when the pump- 
ing contour encircles the resonance line. In the presence of 
a superconducting lead, we find that the pumped charge is 
quantized with the value 2e. 
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Physics of parametria j-alectron pump has attracted 
great attention recentljulla. A classical example of elec- 
tron pump is the Thouless pump facilitated by a traveling 
wave potentialEj. The pumped charge is quantizedEj apd 
can be used as a quantum standard for electric charge.Ej. 
The quantization of pumped charge has also-been stud- 
ied for a large, almost open quantum. dotEj'llj and a 
small, strongly pinched quantum dottZI. In the latter 
case, there exists a resonance line along which the trans- 
mission through the quantum dot is at resonance. The 
pumped charge is quantized if the pumping contour in 
parameterj-space is properly chosen to encircle the reso- 
nance linalll. Recently, we have studied the parametric 
pumping in presence of a superconducting leadt3. At the 
normal conductor-superconductor (NS) interface, an in- 
coming electron-like excitation can be Andreev reflected 
as a hole-like |excitationli3. In contrast to the current 
doubling effectcEl, we found that due to the quantum in- 
terference of direct reflection and the multiple Andreev 
reflection, the pumped current is four times of the value 
when the leads are normal in the weak pumping regime. 
In this paper, we explore the effect of superconducting 
lead on electron pumping in the opposite limit, i.e., we 
study the pumped charge during the pumping cycle in 
the the strong pumping regime. Here the pumped charge 
equal to the pumped current multiplied by the period 
of pumping cycle. Similar to the Ref. |l^, we examine 
the behavior of pumped charge near the resonance line. 
We find that the pumped charge in one pumping cycle 
is quantized with the vaLue of Q = 2e when one of the 
leads is superconducting.Eil 

We consider a parametric pump which consists of a 
double barrier tunneling structure attached to a normal 
left lead and a superconducting right lead. Due to the 
cyclic variation of external parameters xi and X2, the 
adiabatic charge transfer in the presence of a supercon- 
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where r is the period of ( 
dA^L / dx is the injectivitj 
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where the first term is the injectivity p£|electron due to 
the variation of the external parameternjEa, i.e. the par- 
tial density of states (DOS) for an electron coming from 
left lead and exiting the system as an electron, and the 
second term is the injectivity of a hole, i.e. the DOS for 
a hole coming from left lead and exiting the system as an 
electron. Using the Green's theorem, the pumped charge 
can be expressed as surface integral over area A enclosed 
by the path {xi{t),X2{t)) in the parameter spacetl 
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with 



W^^{xi,X2) = Im[ 
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Note that the area A is a, measure of variation of pump- 
ing parameters xi and X2. A is very small in the weak 
pumping limit while remains finite in the strong pumping 
regime. 

For the NS structures, the scattering matrix is de- 
scribed by 2 X 2 matrix S when the Fermi energy is within 
the superconducting gap A. 
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where See (or She) is the scattering amplitude of the in- 
cident electron reflectecLas an electpsjcuior a hole) . Using 
Andreev approximationlij, we havec3cZl 



S = Sii + 5i2(l - RiS22r^RiS2i 



(6) 



where Sp^{E) (/?, 7 = 1,2) is a diagonal 2x2 scatter- 
ing matrix for the double barrier structure with matrix 
element Sp^{E) and S*p^{—E). For instance, we have 
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In Eq.(||) Ri = aa^ is the 2x2 scattering matrix at 
NS interface due to the Andreev reflection with off diag- 
onal matrix element a. Here a = {E — wa/A^ — E^)/ A 
with 1/ = 1 when E > —A and i> — —I when 
E < —A. In Eq.(H), the energy E is measured rela- 
tive to the chemical potential fi of the supierconducting 
lead. Eq.(|^) has a clear physical meanin^Zl. The first 
term is the direct reflection from the normal scatter- 
ing structure arid the second term can be expanded as 
S12R1S21 + S12R1S22R1S21 + ■■■ which is clearly the sum 
of the multiple Andreev reflection in the hybrid struc- 
ture. It is the quantum interference of these two terms 
which gives rise the enhancement of pujBped current in 
the weak pumping regime for NS systenJla. From Eq.(H) 
we obtain the well iaiown expressions for the scattering 
matrix See and Sh(SB 

See{E) = Sn{E) + a'Si2{E)S*22{-E)MeS2i{E) (8) 



pole of the scattering matrixE-0 which works well in one 
dimension or by calculating the dwell time of the mr 
cident electron for two or three dimensional systemsES. 
In the case of double 5 barriers-structure, the energy 
of quasi-bound state is given byll3 E = Er + AE with 
AE ~ ~{kr/a){xi + X2) where Er — — (n-Kja)^ is en- 
ergy of the bound state when the system is isolated. This 
defines a resonance line xi -f 2:2 = — ^ in parameter space 
{x\tX2) along which the transmission is at resonancdiZI. 
Here (5 < is the detuning of the Fermi energy from the 
bound state. 

To show the quantization of charge transfer in the NS 
system, it is useful to recall the calculation of the normal 
case and make the compadscn. In the normal case the 
charge transfer is given byaO 

^— [ dxidX2U^{xi,X2) (15) 
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and 



SheiE) = aSl^{-E)MeS2i{E) 



Me = [l-a^S22{E)S*22{-E)\ 
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The double barrier structure which we consider is mod- 
eled by potential U{y) = Vi5{y + a/2) + V25{y - a/2) 
where Vi and V2 are barrier heights which varys in a 
cyclic fashion to allow the charge pumping. For this sys- 
tem the retardpi Green's function G^{y,y') can be cal- 
culated exactlyEJ. This is done by applying the Dyson's 
equation regarding that any one of the 5-barrier is just a 
perturbation of the remaining system. This way G^{y, y') 
is obtained by applying Dyson's equation twice starting 
from the Green's function of the one-dimensional free 
space. With G^{y,y') we can calculatcj-scattering ma- 
trix exactly from the Fisher-Lee relatione^ 



5/3^ = —5/3^ + ivG^i^^ 



(11) 



where G^^ = G'^{yp, y^) and v — 2k is the electron veloQ- 
ity in the normal lead. For normal structure, we havdlZI 



Sn = [l--iX2~{l + ixi)a^]/D 
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and 



[1 — ix^ — {l+iX2)<J^]/D 



S12 = S21 = XiX2(t/D 



(13) 
(14) 

where D — —(1 — ixi){l — 1x2) + cr^, 2:1^2 — '2kVi^2, 
and a = exp(i/ca). For the double barrier structure, 
the resonant tunneling is mediated by the quasi-bound 
state. When the energy of the incident electron is in line 
with the energy of the quasi-bound state the transmis- 
sion coefficient reaches maximum. The energy of quasi- 
bound states can be determined either by looking at the 



H (-i--2)=/m[^^-K^^] (16) 

The pumped charge in this case has been calculated in 
Ref. n^. In the weak pumping limit, it is easy to show 
that only dxSn contributes to the pumped charge. In 
the strong pumping regime, we will show in the following 
that the contribution from dxSi2 to the pumped charge 
in normal structure is zero. As discussed in detail in 
Ref. |l^, we neglect the smooth energy dependence of xi 
and X2. From Eq.(p^, we obtain the contribution due to 

dxSl2 



Xl^{xi,X2) = Fi{xi,X2)/F2(xi,X2) 



with 



Fi{xi,X2) = -20:12:2(0:1 - 2:2) sin^((5/2) 



(17) 
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^2(^1,2:2) = 2;i.X2 + [xi + 2:2)^ + 2(2:1 -f 2:2) sin (5 

+ 2(1 -2:12:2) (1 -cos (5) (19) 

To compute the surface integral of in Eq.(p7|), it is 
convenient to change the variables from 2:1^2 to p and z: 



and 
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with < p < 00 and — 1 < z < 1. Substituting Eqs.(2C) 
and 
and 



21) into Eqs.(Q and (|l|) and expanding Eqs.(|T§) 
1£) in terms of small 6, we have 



Fi = z(l -z2)^5p3, 



F2^d\l-pf + S'gip,z) 



(22) 
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where q(p, z) (an even function of z) is given in Eq.(8) of 
Ref. 1,1. Since Fi is an odd function of z, the contribution 
due to dxSi2 to the pumped charge is zero. 

Now we follow the same procedure to calculate the 
pumped charge for the NS system. For the parametric 
pumping at zero temperature, we only need the scatter- 
ing matrix at the Fermi level, i.e., at = 0. From Eqs.(H) 
and ([l0|), we see that She is a real quantity and hence 
makes no contribution to the pumped charge in Eq.(||). 
It is straightforward to calculate H^'^ using Eq. (|§) , from 
which we obtain, n^'^(a:i, a:;2) = F^{xi^X2) / Ff{xi,X2) 
where F3 = 4:x\x\{2 — 2 cos 5 -I- a;2sin^) and F4 = 
x\x2 + 2{xi+X2)'^ +^(xi+X2) sin5-|-4(l — a;ia;2)(l — cos(5). 

In Fig.l we plot both 11^'^ and 11^ as well as their 
cross-sections along and perpendicular to the resonance 
line. We see that 11^"^ and 11^ are peaked around the res- 
onance line. Two features are worth noticing. First of all, 
the peak of 11^'^ is much sharper than that of 11^ . This is 
understandable and is due to the resonance nature of NS 
structures near the resonance line. In the Breit-Wigner 
form, the transmission coefficients for normal and NS 
structures-are, respectively |S'2iP = rir2/[(£' — + 
rV4] andEa \She? = ATlTll[i{E-Erf+T\+Tlf, where 
Er is the resonant level, Fi and r2 are the decay widths 
into the left and right lead. Hence |iSheP decays much 
faster away from Er than |S'2ip. The scatteririg ma- 
trix 5*21 and She will appear, respectively, in Eqs.(|^) and 
( |l5| ) implicitly as can be seen from Fisher-Lee relatkai 
Eq.(|n]) and the Dyson equation dx^Gl^ = G^aG^iEil. 
Secondly, the peak height of 11^'^ is four times larger 
than that of 11^. This is precisely due to the construc- 
tive interference of direct reflection and multiple Andreev 
reflectionlia. Now the physics of pumping at resonance is 
clear. For the resonance pumping in the weak pump- 
ing regime, we are looking at the small neighborhood of 
the peak. The area of the neighborhood has to be small 
since it is the weak pumping. The neighborhood has to 
be around the peak with xi ~ X2 since only around the 
peak the transmission coefficient is approximately one. 
As a result, we obtain immediately the pumped charge 
or pumped current of NS structure near the resonance is 
four times of that of corresponding normal structure. In 
the other extreme, for strong pumping, we take a large 
contour enclosing entire resonance line. Since de- 
creases much faster than 11^ away from the peak, it is 
understandable that the pumped charges (the integral of 
n over the area enclosed by the contour) for both nor- 
mal and NS structures are equal, which will be shown 
analytically below. 

After the expansion in powers of 5 in Eqs.(|l^) and ([l9| ) 
and keep the leading orders of 5, we have 



F3 = /[2 +p(-l + z)]{-l + zf{l + zf- (24) 
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lp^(-l + z^) + lp\^l + zy]S' 



So Eq.(||) becomes. 
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^he fact that lim5_>o S^/{x^ + 6^) 
1(m becomes 



(3/8)7r(5(x) 



Q^^^3V-2el\dz^l^-m±£^2e 



■(l-F6z2+z4)5/2 
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Hence the pumped charge for NS system is quantized at 
the same value of that of the normal structure. 

Now we have a better physical picture for the transport 
properties of the NS structure. For the conductance or 
the I-V curve, we need 6*21 or She- For normal structure, 
the current is given by = 2e/h J dE[f{E — eVi) — 
f{E — eV2)]|S'i2p and hence at resonance and at zero 
temperature I^/{Vi - V2) = 2e^/h. For NS struc- 

ture, wehaveBB /^■^ ^ 2e/h J dE{f{E - eVi) ~ f{E + 
eVi))\She\^ and at resonance G^^ = J^^/Vi = ie^/h 
which is the well known doubling of the conductance. For 
pumped charge or pumped current at resonance, how- 
ever, it depends only on dxiSu or dxiSee {E — Q \s 
assumed). Because of the constructive interference be- 
tween direct reflection and multiple Andreev reflection in 
the weak pumping regime, the charge transfer increases 
by a factor of four when one of the lead becomes super- 
conducting. In the strong pumping regime, however, the 
charge transfer is quantized at the value equal to that 
of normal structure, if the pumping contour is chosen 
such that the resonance line is enclosed. The physics be- 
hind this can be understood as follows. In the normal 
case, the contour enclosing the resonance line in the pa- 
rameter space passes through the resonance line at two 
points {xi,X2) = (0, — (5) when the left contact is almost 
closed and (—(5, 0) when the right contact is almost closed. 
When passing through those two points, the resonance 
level of the dot crosses the Fermi energy. At each cross- 
ing, the occupation of the level changes, and two electrons 
with opposite spin enter or exit the region between the 
barriers. Since one of the tunnel barriers has zero conduc- 
tance at those points, it is clear that the electrons must 
have tunneled through the other contact upon entering 
or leaving the quantum dot. Hence, in the pumping cy- 
cle, electrons are shuttled pairwise through the dot. In 
the presence of superconducting lead, the resonance level 
(both the energy and the width) is exactly the same as 
that of normal case since the scattering matrix is given 
by She = i|5'i2|V(l + \S22?) when E = Q. Therefore the 
same argument applies to the superconducting case and 
the quantization unit is 2e. Note that our statement is 
only valid when the electron interaction is neglected. For 
the case of two normal-metal contacts, if interactions are 
included the quantization will remain, but now the quan- 
tum is only e: Only one electron at a time can enter the 
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region between the barriers; addition of a second electron 
is forbidden by Coulomb blockade. In the presence of su- 
perconducting lead, since the Andreev reflection requires 
two electrons with opposite spin in order to produce the 
supercurrent, it seems that the pumping is not allowed in 
the strong pumping regime due to Coulomb blockade. In 
this paper, we have also neglected the effect of tempera- 
ture and the effect of inelastic scattering. As discussed in 
Ref. |l^ the temperature will destroy the quantization of 
the pumped charge. When inelastic channel is present an 
additional physicaLjpcchanism for an incoherent pump 
effect will show upE^. 
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FIG. 1. (a). The integrand U of Eqs.(3) and (15) as a 
function of xi and X2 for 5 = —0.2. For illustrating purpose, 
the origin of n^(a;i,a;2) has been shifted by (0.1,0.1). (b). 
The cross-section of II along the resonance line xi+X2 ~ —S. 
Solid line: R'^^; dotted line: n'^ . Inset: the cross-section of 11 
along the direction xi — X2 = co which is perpendicular to the 
resonance line. Left inset: co = 0.01; right inset: cq — —0.042. 
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